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ABSTRACT 

Recent work by Risaliti et al. (1999) suggests that more than half of all Seyfert 2 
galaxies in the local universe are Compton-thick (A'h > lO^'^cm"^ ). This has impli- 
cations for AGN synthesis models for the X-ray background (XRB), the flexibility 
of which for the inclusion of large numbers of high-z type 2 sources we examine 
here. We highlight the importance of Compton down-scattering in determining the 
individual source spectra and the fit to the XRB spectrum, and demonstrate how 
parameter space 'opens up' considerably if a super-solar iron abundance is assumed 
for the absorbing material. This is illustrated with a model which satisfies the present 
constraints, but which predicts substantial numbers of type 2 sources at the faint 
flux levels soon to be probed for the first time by the Chandra and XMM missions. 
We demonstrate also how a strong negative K-correction facilitates the detection of 
sources with 10~^'* < A^h < 10^^ cm~^ out to the highest redshifts at which they could 
plausibly exist. 
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1 INTRODUCTION 

There have been numerous attempts to account for the ori- 
gin of the cosmic X-ray background (XRB) with populations 
of active galactic nuclei (AGN) integrated over cosmic time 
(Setti & Woltjer 1989; Madau et al. 1994; Celotti et al. 1995; 
Comastri et al. 1995). In this case, Fabian & Iwasawa (1999) 
have recently shown that the XRB sources plausibly account 
for all of the accretion power in the Universe, and that most 
of it is intrinsically absorbed. The XRB is thus an impor- 
tant probe of the evolution of the black hole population with 
cosmic epoch. It is of interest to know how this relates to 
the evolution of star formation in galaxies, which work at 
sub-mm wavelengths is beginning to constrain (eg. Hughes 
et al. 1998). 

Within the framework of unified models for AGN, XRB 
synthesis models require as inputs the X-ray spectra of the 
unabsorbed (type 1) and absorbed (type 2) sources, and 
the local X-ray luminosity function (XLF) for a particular 
class of source together with its cosmological evolution. It 
has been common practice to adopt an XLF for the type 1 
sources, and to convolve it with an Ah distribution prescrib- 
ing the relative numbers of sources with different amounts 
of absorption. The form of the Ah distribution is, however, 
the major uncertainty in current models; the earliest works 
assumed it to be independent of both luminosity and red- 
shift, and consistent with the poorly-constrained local esti- 



mates of the Seyfert 2/Seyfert 1 number ratio. The latter 
lies between ~1.6-3 according to Huchra & Burg (1992) 
for an optically-selected Seyfert sample, but other studies 
constrain it to lie between 2 and 10 (Osterbrock & Shaw 
1988; Goodrich, Veilleux & HiU 1994). More recently, there 
have been attempts to determine the Ah distribution di- 
rectly: from a large sample of Seyfert 2 galaxies, Bassani 
et al. (1998) found a mean log Ah of 23.5 with 23-30 per 
cent of sources having log Ah > 24. Their distribution is, 
however, biased in favour of low Ah values, as they include 
all Seyfert 2s for which a 'reliable' hard X-ray spectrum is 
available in the literature. Risaliti et al. (1999) overcame 
this and other selection biases by using the [OIIIJA5007 flux 
to select a sub-sample of Seyfert 2 galaxies from the sample 
of Maiolino & Rieke (1995); they found an Ah distribution 
even more strongly dominated by heavily absorbed sources, 
with 75 per cent having log Ah > 23 and about half having 
log A^'h > 24. They found no correlation between Ah and 
the intrinsic luminosity of the nuclear source. 

Attempts to incorporate some of this new information 
into the synthesis models have recently been made (Gilli et 
al. 1999a). We emphasize that the above Ah distributions 
are determined locally, and that the situation at higher red- 
shifts may be quite different. To this end, we examine in 
the present work how much flexibility there is within the 
models for accommodating type 2 sources in substantially 
larger numbers at high redshift, at flux levels within reach of 
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the forthcoming generation of X-ray observatories (Chandra 
and XMM). In so doing, we draw attention to the role played 
by Compton down-scattering in determining the spectra of 
the individual sources and hence the quality of the fit to 
the XRB spectrum. Any model which includes large num- 
bers of Compton-thick sources must clearly incorporate this 
feature; it was omitted by Comastri et al. (1995), and (as 
pointed out by Celotti et al. 1995) handled erroneously by 
Madau et al. (1994). Only Celotti et al. (1995) appear to 
have treated it properly, and they found a very restricted 
range of parameter space for acceptable models. We demon- 
strate that the quality of the fit to the XRB spectrum can 
be improved, and the available parameter space enlarged, by 
positing a super-solar metallicity for the absorbing material. 

We adopt the cosmological parameters Ho ~ 
50kms~^ Mpc"'^ and go = 0.5 throughout. 



2 THE MODEL 

2.1 The AGN XLF and evolution 

In common with Madau et al. (1994) and Celotti et al. (1995) 
we adopt the Piccinotti et al. (1982) 2-10 keV XLF for 
local type 1 sources: n(L44) = nol/^/'^^, where L44 is 
the 2-10 keV luminosity in units of 10'*'' erg s~^ and no — 
2.68 X 10"^Mpc-^(10''*ergs-i)-\ between Lx.min = 7 x 
10''^ergs~* and Lx,7nax = WOLx.min- Others (eg. Gilh 
et al. 1999a, 1999b, Comastri et al. 1995), use the 0.3- 
3.5 keV XLF for broad line objects (QSOs and Seyfert Is) 
derived from ROSAT and (Einstein) EMSS data (Jones et 
al. 1997 and references therein). Whilst the ROSAT-I-EMSS 
XLF extends to higher luminosities than that of Piccinotti 
et al. (1982), the integrated number density derived from 
the latter is larger by an order of magnitude, showing 
that there are many 2-10 keV sources absent from the 0.3- 
3.5 keV band. It is unlikely that this discrepancy can be 
wholly ascribed to the contamination of the Piccinotti et al. 
sample by type 2 sources (of which there are 7, compared 
with 18 of type 1). 

We adopt pure luminosity evolution of the form (1 + 
z)^ out to a cut-off redshift Zcut, with no further evolution 
thereafter out to a redshift Zmax- This is consistent with 
the findings of Jones et al. (1997) who derive Zcut ~ 1.4 
and p ~ 3. We note that Miyaji et al. (1998) have recently 
determined the XLF in the observed 0.5-2.0 keV band for 
all AGN (ie. not just the type 1 objects) by combining the 
existing ROSAT survey data. They propose a luminosity- 
dependent density evolution model, but do not provide a 
break-down according to source classification, so the effects 
of intrinsic absorption render their XLF of limited use for 
synthesis models. 



2.2 Source spectra 

For the primary type 1 X-ray spectrum, we follow Madau et 
al. (1994) and assume the form I{E) oc E'^'^ exp{-E / Ec) 
with Ec fixed at 360 keV. To this is added a component re- 
fiected from an accretion disc subtending 2-k sr at the source, 
for a fixed inclination angle of 60°. Type 2 spectra are gen- 
erated by feeding this primary spectrum to a Monte-Carlo 



code which computes the transmission through a homoge- 
neous sphere of cold material with a radius given by the col- 
umn density A'h . Both photoelectric absorption from neutral 
material (of solar composition unless explicitly stated other- 
wise) and Compton down-scattering are modelled, the latter 
being described by the Klein-Nishina cross-section. Investi- 
gations with the obscuring material in the form of shells and 
slabs showed that the output spectra are largely indepen- 
dent of the assumed geometry. Finally, a component equal 
to 2 per cent of the primary incident power is added to the 
type 2 spectra to represent the flux scattered into the line 
of sight by electrons in the warm, ionized medium above 
the axis of the obscuring torus. The Monte-Carlo spectra 
closely match those generated independently by Celotti et 
al. (1995) (A. Celotti, private communication). 

We note that Celotti et al. (1995) assumed the sources 
to be compact and in thermal pair equilibrium, as a result of 
which the cut-off energy Ec evolves with source luminosity. 
For simplicity, we do not model this here as it can be shown 
that if a range of Ec values is employed, the mean must 
be ~ 300 keV in order to match the slope of the observed 
XRB (Gruber et al. 1992) between 30 and 100 keV (which 
is almost entirely intrinsically unabsorbed). 

The A^H distribution is parameterised as A''(logAi'H) oc 
(logiVn)'' from N^.rain = 10^"cm-2 to Nh.* = 
10^^ cm~^ and from A^h,* to Ai'H,max for sources of types 
1 and 2, respectively. In a departure from previous works, 
we investigate the effects of a redshift-dependent /3, but do 
not incorporate any explicit luminosity dependence. The re- 
sults of Risaliti et al. (1999) are consistent with the latter 
assumption (at least locally), and with /3 ~ 8. 

In Fig. N we show the spectra generated by our Monte- 
Carlo code for various A^h values, assuming solar metallic- 
ity for the absorbing material. It is clear that all the spec- 
tra peak in EI(E) at E < 30keV. Thus, regardless of the 
details of the Ah distribution and of the form of the cos- 
mological evolution, there is little scope for generating an 
XRB spectrum peak at the observed energy of 30 keV, espe- 
cially if large numbers of high redshift sources are included. 
Within the formalism of the present synthesis model, the 
only means of 'opening up' the parameter space is to posit 
a super-solar metallicity for the absorbing material. Specif- 
ically, by raising the iron abundance the photoelectric ab- 
sorption can be increased whilst not substantially affecting 
the amount of Compton down-scattering, thereby shifting 
the peaks in the EI{E) spectra to higher energies (the po- 
sitions of the peaks are relatively insensitive to the value of 
Ec). This is illustrated in Fig. H for an iron abundance of 
5 X solar. If Compton down-scattering is 'turned off' so that 
photoelectric absorption remains the only source of opacity, 
this particular difficulty does not arise, as Fig. M shows. 



3 RESULTS 

We discuss here the model fits to both the XRB spectrum 
and the source counts (log N — log S) in the several bands, 
for various values of the evolutionary parameters Zcut, Zmax 
and p, and different A^h distributions. We parameterise the 
latter with the following redshift-dependence for /3: 
/3 = /3i for < z < 0.5 
l3 = (3-2 for 0.5 < z < 1.0 
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Figure 1. For the primary input spectrum described 
in the text (dashed line), the solid lines are the type 
2 spectra produced by the Monte-Carlo code (plus 
the 2 per cent component scattered by warm ion- 
ized material) for the column densities (from left to 

right) TVh = 10^1.25^^022.25 ;^023.25^;^024.25^ ^024.75 

10^^ '^^ cm"^ . Note that all the spectra peak &t E < 30keV, 
as the position of dotted vertical line makes clear. 
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Figure 2. As in Fig. |l|, but for an iron abundance of 5x solar. 
Note how the peaks in the spectra have been shifted to > 
30keV, which is crucial for reproducing the observed peak 
in the XRB spectrum at 30 keV. 

/3 = /33 for 1.0 < 2 < 2.0 
/? = /34 for 2.0 < z < 3.0 

l3 = for 3.0 < 2 < Z„^ax 

We take A'^e.max ~ 10^^ cm~^ , as Figs. ^ and ^ show that 
there is essentially no direct flux (and hence little contri- 
bution to the XRB) from more heavily obscured sources 
(but Fig. ^ shows that this is not so when Compton down- 
scattering is neglected). Acceptable fits to the observed 
XRB spectrum as parameterised above 3keV by Gruber et 
al. (1992) were defined as satisfying the criteria given by 
Celotti et al. (1995). As stated therein, the model spectrum 
is first re-normalized to the observed XRB at 30 keV by be- 
tween ~ 0.7 and 1.5, to take into account uncertainty in the 
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Figure 3. As in Fig. |l| but with Compton down-scattering 
'turned off; note how the spectral peaks of the Compton- 
thick sources are shifted well above 30keV, thus creating 
extra parameter space for XRB synthesis models which ne- 
glect this important effect. 

normalization of the type 1 XLF, as set by the assumed value 
of A'h,* and the level type 2 contamination in the Piccinotti 
et al. (1982) XLF. 

3.1 Models with solar iron abundance 

With solar iron abundance and /3i_5 = 8, acceptable models 
can be found for Zcut ~1.2"1.4, p ~2.5-2.9 and Zmax —3.5- 
4.0 (ie. within the 68 per cent confidence limits derived by 
Jones et al. 1997). There is no freedom for the inclusion 
of larger numbers of type 2 sources at high redshift. The 
predicted local emissivities in the 2-10 keV band lie at or 
just beyond the upper end of the range ~ (3.5 — 4.8) x 
10^**h5o erg s"^ Mpc"^ derived by Miyaji et al. (1994), from 
the cross-correlation of the XRB measured with HEAO-1 
and IRAS galaxies. The models resolve ~ 95 — 100 per 
cent of the 2-10 keV XRB intensity measured by Mar- 
shall et al. (1980), which can be compared with an esti- 
mated contribution from clusters of galaxies of ~ 4.7 per 
cent. We obtained the latter figure by integrating the Ebel- 
ing et al. (1997) cluster XLF over 10^^ < L(2 - lOkeV) < 
10'*® ergs~^ , assuming the David et al. (1993) Lx — T rela- 
tion for a thermal bremsstrahlung spectrum and no cosmo- 
logical evolution out to a maximum redshift z=4. The de- 
rived log — log S relation in the 2-10 keV band is consistent 
with the constraints from Piccinotti et al. (1982), Butcher 
et al. (1997) and with the recent measurements down to 
5(2 - lOkeV) = lO-^^ergcm-^s"^ from the ASCA Large 
Sky Survey (Ueda et al. 1999). 

3.2 Models with super-solar iron abundance 

As discussed in section 2, much more parameter space is 
available to satisfactory models if the metallicity of the ab- 
sorbing material is increased. In particular, there is scope for 
including substantially more type 2 sources at high-redshift. 
This is illustrated in Fig. ^ where we show the percentage 
deviation of the model from the Gruber (1992) fit to the 
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observed spectrum for various TVh distributions and for two 
different iron abundances. We consider models constructed 
using an iron abundance of 5 times solar, which is perhaps 
an extreme value but which should bracket the likely range. 
Rather than provide an extensive tour of parameter space, 
we illustrate what is possible with an arbitrarily chosen 
model for which /3i_3 = 8, /34,5 = 16, Zcut = 1.3, Zmax = 4.0 
and p = 2.5. After renormalization at 30 keV (by a factor 
of 0.93), this model resolves 90 per cent of the Marshall et 
al. (1980) 2-10 keV XRB intensity, with a 2-10 keV emissiv- 
ity of 4.30 X lO^^^hso erg s"^ Mpc"^. [Formally acceptable fits 
with some of the extreme /? values of Fig. ^ can be obtained 
with a solar iron abundance by substantially increasing the 
value of -Ec and renormalizing appropriately; such models, 
however, yield XRB spectra with slopes which are too flat 
above 30keV, as mentioned in section 2.2]. 

Fig. ^ shows the fit to the Gruber et al. (1992) param- 
eterization of XRB spectrum, along with the contributions 
of the different classes of source. Figs. I and show the 
model predicted source counts in the hard (2-10 keV) and 
soft (0.5-2 keV) bands respectively, along with the current 
observational constraints. Fig. ^ saturates at a flux level al- 
most two orders of magnitude lower than the models pro- 
posed by Celotti et al. (1995), in part because the latter au- 
thors used a lower Zmax~'^-Q (and a correspondingly higher 
2ciit=2.0). Compton-thin type 2 sources are dominant in the 
interval lO"^'^ < 5(2 - lOkeV) < 10"^" erg cm'^ s'^ and 
the surface density of Compton-thick sources in 10~^® < 
5(2 - lOkeV) < lO-^'^ergcm-^s"^ is ~ 1600 deg-^ This 
is of interest as the forthcoming Chandra and XMM mis- 
sions will probe down to 2-10 keV flux levels of ~ 10~^® and 
~ 10~^^ erg cm"^ s~^ , respectively. For comparison with 
the preliminary BeppoSAX High Energy Large Area Sur- 
vey (HELLAS) results (Comastri et al. 1999), we show in 
Fig. y the source counts in the 5-10 keV band; the discrep- 
ancy at intermediate fluxes may be due to the contribution 
of non-AGN sources or to residual calibration problems, as 
described by the latter authors. Because of its potential rele- 
vance to future missions, we show in Fig. ^ the source counts 
for the 10-20 keV band. Fig. ^ shows the fractions of sources 
of each type as a functions of 5(10 - 20 keV), 5(2 - 10 keV) 
and 5(0.5- 2.0 keV). 

The discrepancy between the model-predicted and ob- 
served source counts in the 0.5-2.0 keV band perhaps indi- 
cates that our modelling of the soft X-ray spectra of the 
sources is unrealistic. Indeed, it is partly because there are 
large uncertainties in the spectral shape of the soft X-ray 
emission in AGN that we follow Celotti et al. (1995) in ne- 
glecting the soft X-ray band when assessing the quality of 
the fit to the XRB spectrum. There is some scope for adjust- 
ing the model counts at faint 0.5-2.0 keV fluxes in Fig. Mby 
changing the scattering fraction for the type 2 sources from 
the figure of 2 per cent used here. That there might be some 
dispersion in this quantity is suggested by ROSAT HRI ob- 
servations of some hyperluminous IRAS galaxies, where the 
upper limits on the 0.1-2.4 keV fluxes imply scattering frac- 
tions below ~ 0.5 per cent (Fabian et al. 1996; Wilman et 
al. 1998). 
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Figure 4. For iron abundances of solar (triangles) and 5 
times solar (asterisks), each of the panels shows, for a differ- 
ent TVh distribution, the percentage deviation of the model 
from the observed (Gruber 1992) XRB spectrum at energies 
between 10 and 100 keV, after normalization at 30keV. In 
each panel, the models shown provide the best-flt to the 
spectrum for Zcut — 1.3, Zmax — 4.0 and p ~ 2.5. The 
±10 per cent limits for acceptable models and the (3 pa- 
rameters of the A^H distributions are as indicated. There is 
a larger range of permitted A'h distributions with the higher 
iron abundance. 



3.3 Detecting Compton-thick sources at high 
redshift 

The 'kink' which appears at a flux of 10"^'^ ergcm"^ s~^ in 
the contribution of Compton-thick type 2 sources to Fig. ^ is 
due to a negative K-correction which applies to the spectra 
of these sources (a similar feature ispresent in the Compton- 
thin type 2 contribution to Fig. Indeed, Fig. ^ shows 
that the spectrum of a source with A''h ~ 10^*'^ cm"'^ rises 
very steeply above lOkeV, and the redshifting of this rest- 
frame 30 keV peak down into the observed 2-10 keV band 
can thus more than offset the D^'^ cosmological dimming. 
The extent of the effect is illustrated in Fig. ^ where 
we show the observed flux from such a source as a func- 
tion of redshift, for several values of the scattering frac- 
tion. It is thus apparent that deep Chandra and XMM 
pointings will be able to detect the high luminosity end 
[L(2 — lOkeV) > 10*^ ergs"'^ in the unabsorbed rest-frame] 
of the XLF for Nn < 10^^ cm'^ (above which A^h the K- 
correction is much smaller) out to the highest redshifts at 
which such sources could plausibly exist. This is analogous to 
the way in which the sub-mm-detectability of star-forming 
galaxies with Arp 220-like rest-frame far-infrared SEDs is 
greatly enhanced between redshifts 2-5 (Blain & Longair 
1993). 

We note that the models are much less sensitive to the 
value of Zmax than they are to either p or Zcut- The existence 
of sources at z » Zmax thus cannot be ruled out. 



© 0000 RAS, MNRAS 000, 000-000 



Fitting the spectrum of the X-ray background: the effects of high metallicity absorption 5 



o 




° I 10 100 



E (keV) 

Figure 5. Comparison of the XRB spectrum of the model de- 
scribed in the text for an iron abundance of 5 x solar (dashed 
line) with the analytical fit of Gruber et al. (1992) to the ob- 
served spectrum (solid line). The contributions to the model 
from type 1 sources (dot-dashed line), and the Compton-thin 
(dotted line) and Compton-thick (three dots-dashed line) 
type 2 sources are as shown. 




Figure 6. Predicted log A*' - log S in the 2-10 keV band of 
the model described in section 3.2. The solid line shows the 
sum of all the source counts, the dashed line the type 1 
sources with 20 < log A'h < 22, the dash-dot line type 2 
sources with 22 < log A^h < 24 and the dotted line those 
sources with 24 < logiVn < 25. The cross is the HEAO-1 
A2 AGN surface density from Piccinotti et al. (1982), the 
asterisk the result of GINGA fluctuation analysis after the 
removal of the cluster contribution (Butcher et al. 1997) and 
the filled circles represent the results of the ASCA Large Sky 
Survey (Ueda et al. 1999). 

4 SUMMARY AND CONCLUSIONS 

We have re-examined AGN synthesis models for the XRB 
in the light of evidence suggesting that more than half of all 
Seyfert 2 galaxies in the local universe are Compton-thick 
(Risaliti et al. 1999). Compton down-scattering and the iron 
abundance of the absorbing material influence the appear- 
ance of the source spectra around 30 keV in ways which sig- 




Figure 7. As in Fig. ^ but for the 0.5-2.0 keV band. The open 
circles are from the ultradeep (1 Msec) ROSAT HRI obser- 
vation in the Lockman Hole (Hasinger et al. 1998), and the 
solid points are from shallower ROSAT surveys (Hasinger et 
al. 1993). 
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Figure 8. As in Fig. |6| but for the 5-10 keV band. The points 
are preliminary BeppoSAX HELLAS results (Comastri et 
al. 1999). 

nificantly affect the fit to the XRB spectrum, as summarised 
below. 

(1) It is possible to fit the XRB spectrum assuming a 
solar iron abundance, but there is no scope for accommo- 
dating an increased type 2/type 1 ratio at higher redshifts. 

(2) With an iron abundance of 5x solar, the fit to the 
XRB spectrum is better and much more parameter space 
becomes available. We illustrated this with an arbitrarily 
chosen model in which the assumed type 2/type 1 ratio 
increases from a local value of 3.8 to 9.7 above z = 2. 
This model predicts that the forthcoming X-ray missions 
will uncover type 2 sources in substantial numbers below 
S{2 — lOkeV) = lO"^** erg cm"'^ s~^ . Furthermore, a strong 
negative K-correction facilitates the detection of Compton- 
thick sources with lO^"* < A''h < lO^'^ cm"^ out to high 
redshift. 

As our aim was to demonstrate what is possible within 
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Figure 9. As in Fig. ^ but for the 10-20 keV band, potentially 
important for future missions. 
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Figure 10. The fraction of type 1 (solid line), Compton- 
thin (dashed line) and Compton-thick (dot-dashed line) 
type 2 sources, are shown as functions of 5(10 — 20keV) 
5'(2 - lOkeV) and 5(0.5 - 2.0 keV) for the model described 
in section 3.2, with an Fe abundance of 5x solar. 



the models, we do not discuss the cosmological constraints 
bearing upon the existence of large iron abundances in high- 
redshift AGN. We note though that the strength of the 
NVA1240 emission-line in some quasars demands metallic- 
ities of at least 5Zq (Ferland et al. 1996), consistent with 
models of early-epoch rapid star formation. In a recent re- 
view, Hamann and Ferland (1999) concluded that there is 



Figure 11. For a completely unabsorbed (rest-frame) 2- 
lOkeV luminosity of 10*^ ergs~^ , the curves show the ob- 
served 2-10 keV flux as a function of redshift from a source 
with A^H ~ lO'^'*'^ cm~'^ , for scattering fractions of (from 
top to bottom) 5, 1, 0.1 and 0.01 per cent, respectively. The 
strong negative K-correction acting upon the 30keV rest- 
frame peak is clearly evident. 



a growing consensus for super-solar metallicities in quasars 
out to 2: > 4 (typically at the level of one to a few times 
solar), with the suggestion of higher metallicities in more 
luminous quasars and with no evidence for a decline at the 
highest redshifts. Concerning iron itself, they noted that a 
super-solar Fe/Mg abundance out to 2: > 4 is one possible 
solution to an energy budget problem with the broad Fell 
emission lines. We also note that studies of the fluorescent 
iron Ka line in some Seyfert 1 galaxies [see eg. Reynolds 
et al. (1995) and Lee et al. (1999)] are consistent with iron 
abundances of ~ 2x solar. 

A super-solar metallicity can, however, plausibly be 
accommodated within the absorption model of Fabian et 
al. (1998) for the sources which constitute the XRB. The 
latter authors noted the need for an obscuring medium with 
a large covering fraction as seen from the central source, in 
order to generate the high fraction of obscured AGN which 
early XRB models recognised as being necessary (Setti & 
Woltjer 1989). They thus proposed that nuclear starbursts 
obscure the XRB, with the continual injection of energy from 
supernovae serving to keep the medium well 'inflated'. With- 
out this mechanism it would rapidly dissipate into the 'stan- 
dard' dusty torus invoked by unifled schemes, leading to too 
low a covering fraction. These same supernovae could supply 
and subsequently maintain (by, for example, preventing de- 
pletion onto dust grains) a high gas phase iron abundance. A 
high metallicity would also reduce the mass of gas required 
for a given amount of absorption. 
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